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Abstract 



The relativistic mean field (RMF) model is applied to a system of nucle- 
ons and a K meson, interacting via scalar and vector boson fields. The model 
incorporates the standard RMF phenomenology for bound nucleons and, for 
the K meson, it relates to low-energy and K~ atom phenomenology. 
Deeply bound K nuclear states are generated dynamically across the periodic 
table and are exhibited for 12 C and 16 over a wide range of binding energies. 
Substantial polarization of the core nucleus is found for these light nuclei. 
Absorption modes are also included dynamically, considering explicitly both 
the resulting compressed nuclear density and the reduced phase space for K 
absorption from deeply bound states. The behavior of the calculated width 
as function of the K binding energy is studied in order to explore limits on 
the possible existence of narrow K nuclear states. 
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I. INTRODUCTION 



The KN interaction near threshold is strongly attractive, in agreement with the existence 
of the unstable bound state A(1405) below the K~p threshold. The ^-nucleus interaction is 
also strongly attractive, as derived from the strong-interaction shifts and widths in kaonic- 
atom levels across the periodic table [1] . It is not established yet how strong the ^-nucleus 
potential is: 'deep' (150-200 MeV [2,3]) or relatively 'shallow' (50-60 MeV [4,5])? Is it 
possible to bind strongly K mesons in nuclei and are such potentially deep bound states 
sufficiently narrow to allow observation and identification? These issues have received con- 
siderable phenomenological and theoretical attention recently [4-8], and some experimental 
evidence for candidate states in the (K~ top , n) and {K^ top , p) reactions on 4 He (KEK-PS E471, 
[9,10] respectively) and in the (K~,n) in- flight reaction on 16 (BNL-AGS, parasite E930 
[11]) has been presented very recently. New experiments have been approved at KEK, using 
(K~,N) reactions to search for K nuclear bound states [12]. (K~,ir~) reactions [13] were 
also suggested in this context. 

A prime concern in searching for K nuclear bound states is the anticipated large width 
due to pionic conversion modes on a single nucleon: 



with thresholds about 100 MeV and 180 MeV, respectively, below the KN total mass, and 
also due to non-pionic multi-nucleon absorption modes, say 



with thresholds about m w = 140 MeV lower than the single-nucleon thresholds. The branch- 
ing ratios in parentheses are known from bubble-chamber experiments [14]. 

The aim of the present work is to study dynamical effects for K nuclear states in the 
range of binding energy ~ 100 — 200 MeV [9-11] and in particular the width anticipated 
for such deeply bound states. The relatively shallow chirally-motivated ^-nucleus potentials 
[4,5] which followed the microscopic construction by Ramos and Oset [15] are of no use in 
this context, since they cannot yield binding energy greater than the potential depth of 
about 50 MeV. One must therefore depart from the microscopic approach in favor of a more 
phenomenologically inclined model which is constrained by data other than two-body KN 
observables. The theoretical framework here adopted is the relativistic mean field (RMF) 
model for a system of nucleons and one K meson interacting through the exchange of scalar 
(er) and vector (cu) boson fields which are treated in the mean-field approximation. By 
allowing the K to polarize the nucleons, and vice versa, this dynamical calculation is made 
self consistent. K absorption modes are included within a tp optical-model approach, where 
the density p plays a dynamical role, and the constant t which is constrained near threshold 
by fT~-atom data follows the phase-space reduction in reactions Eqs. (1,2) for a deeply 
bound K. A wide range of binding energies may be explored in the RMF calculation simply 
by scanning over the coupling constants of the K meson to the a and uo boson fields. Detailed 
calculations were done by us across the periodic table. In this Letter we demonstrate the 
essential points and conclusions for 12 C and 16 where the dynamical polarization effects 
are extremely important for the energies of the K bound states as well as for their widths. 



KN^irZ, ttA (-80%) , 



(1) 



KNN — >• YN (-20%) , 



(2) 
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Akaishi et al. [7,8] too have found large polarization effects in lighter nuclei using few-body 
variational techniques. The RMF is a systematic approach used across the periodic table 
beyond the very light elements explored by other techniques, and it can be used also to 
study multi--K" configurations and to explore the K condensation limit [16,17]. Similar RMF 
calculations have been recently reported for N states in nuclei [18,19]. 



II. METHODOLOGY 

In the calculations described below, the standard RMF Lagrangian C^- with the linear (L) 
parameterization of Horowitz and Serot [20] as well as the nonlinear (NL) parameterization 
due to Sharma et al. [21] are used for the description of the nucleonic sector. The (anti)kaon 
interaction with the nuclear medium is incorporated by adding to C^- the Lagrangian density 
Ck. [16,17]: 

C K = V*KVi*K - m\KK - g uK m K oKK . (3) 

The covariant derivative = <9 M + ig^K^n describes the coupling of the (anti)kaon to the 
vector meson u. The vector field u is then associated with a conserved current. The coupling 
of the (anti)kaon to the isovector p meson is here excluded due to considering N = Z nuclear 
cores in this initial report. 

Whereas adding C K to the original Lagrangian C_m does not affect the form of the 
corresponding Dirac equation for nucleons, the presence of K leads to additional source 
terms in the equations of motion for the meson fields a and uj q to which the K couples: 

(- A + mty a = -g aN p s - g aK m K KK + (-g 2 a 2 - g 3 cr 3 ) , (4) 



(-A + m£) coo = +guNpv - ^9^k{^k + g U KW Q )KK , (5) 

where ojk is the K energy in the nuclear medium: 

= \Jm 2 K + g aK m K o + p 2 K - g uK Wo , 

and ps and py denote the nuclear scalar and vector densities, respectively. Adding K to 
the nuclear system affects the scalar and vector potentials which enter the Dirac equation 
for nucleons. This leads to the rearrangement, or polarization of the nuclear core in the 
presence of K. 

In order to preserve the connection to previous studies of kaonic atoms, the Klein Gordon 
(KG) equation of motion for the K is written in the form [5]: 

A - 2p(B + V opt + V c ) + (V c + Bf] K = (h = c=l). (6) 

Here, V c denotes the static Coulomb potential for the K~ , p is the K- nucleus reduced mass 
and B = Bk + ^k/1 is the complex binding energy. The real part of the K optical potential 
V ov t is then given by 
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n eK ^ J ( ls-v-^-) , (7) 

where S = g a K& and V = g^K^o are the scalar and vector potentials due to the a and uj 
mean fields, respectively. 

Since the RMF approach does not address the imaginary part of the potential, Im V op t 
was taken in a phenomenological tp form, where its depth was fitted to the K~ atomic 
data [3]. Note that p in the present calculations is no longer a static nuclear density, but 
is a dynamical entity affected by the K interacting with the nucleons via boson fields. The 
resulting compressed nuclear density leads to increased widths, particularly for deeply bound 
states. On the other hand, the phase space available for the decay products is reduced for 
deeply bound states, which will act to decrease the calculated widths. Thus, suppression 
factors multiplying Im V opt were introduced from phase-space considerations, taking into 
account the binding energy of the kaon for the initial decaying state, and assuming two-body 
final-state kinematics for the decay products. Two absorption channels were considered. In 
the first, Eq. (1), a KN initial state decays into a nY final state. The corresponding 
density-independent suppression factor is given by 



f _ M" 01 
Ji 



[Ml - (m, + m Y ) 2 ][Mf - (my - m„) 2 ] _ ^ _ ^ 



1 Mf \| [Ml - (m, + my) 2 ] [M 2 , - (my - m„) 2 ] 

where M i = m K + mjy, M\ = M m — Bk- In the second absorption channel, Eq. (2), 
a KNN initial state decays into a YN final state. The corresponding suppression factor is 
given by 



f _ M 02 

H 



[M| - (m N + mYnM ] - (m Y - m N f] ^ ^ ^ (g) 



^1 ^ [M 2 2 - (mjv + my) 2 ] [M 2 2 - (m y - mjv) 

where M 02 = + 2m N , M 2 = M 02 —B K . Although multi-nucleon absorption modes 
are often modeled to have a power-law p a (a > 1) density dependence, our comprehensive 
i^~-atom fits [2,3] are satisfied with a ~ 1. We therefore assume a — 1 in this exploratory 
work also for this second absorption channel,* which means that / 2 too is independent of 
density. We comment below on the effect of a possible density dependence of / 2 , reflecting 
perhaps a p 2 dependence of the non-pionic decay mode (2) at high densities. 

Since E final states dominate these channels [14] the hyperon Y was here taken as Y = S. 
Allowing A hyperons would foremost add conversion width to K states bound in the region 
Bk ~ 100 — 180 MeV. For the combined suppression factor we assumed a mixture of 80% 
mesonic decay and 20% nonmesonic decay [14], i.e. 

/ = 0.8 h + 0.2 f 2 ■ (10) 

This suppression factor is plotted as function of Bk in the upper part of Fig. 1, where a 
residual value of / = 0.02, when both f\ and f 2 vanish, was assumed. 



*a lucid theoretical discussion of this point, with allowance for a = 1, is due to Koltun [22] 
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The coupled system of equations for nucleons and for the electromagnetic vector field 
A , and for the mean fields a and c^o, Eqs. (4,5) above, as well as the KG equation (6) for 
K~ were solved self-consistently using an iterative procedure. Obviously, the requirement 
of self-consistency is crucial for the proper evaluation of the dynamical effects of the K 
on the nuclear core and vice versa. We note that self-consistency is not imposed here on 
the final-state hadrons which only enter through their on-shell masses used in the phase- 
space suppression factors given above. For the main 7rS decay channel it is likely that the 
attraction provided by the pion within a dynamical calculation [15] is largely cancelled by 
the nuclear repulsion deduced phenomenologically for E hyperons [23-25]. 

III. RESULTS 

The main objective of the present calculations of ^-nucleus bound states was to es- 
tablish correlations between various observables such as the K binding energy, width and 
macroscopic nuclear properties. In particular we aimed at covering a wide range of bind- 
ing energies in order to evaluate widths of possible strongly bound K states. Furthermore, 
in order to study effects of the nuclear polarization, we calculated rms radii and average 
densities of the nuclei involved and, in some cases, also single particle energies. Extensive 
calculations were made for 12 C and 16 0, nuclei that had been discussed earlier in the context 
of strongly bound K states [6,11]. Additionally, more restricted calculations were made for 
40 Ca and 208 Pb. 

In a preliminary test we performed dynamical calculations for K~ atomic Is states, 
which produced only negligibly small polarization effects, thus validating previous analyses of 
kaonic atom data. The empirical values and g^x, as found from a fit to kaonic atom data 
[3], were therefore used as a starting point for calculations. A full dynamical calculation was 
then made for K~ nuclear states starting from the tp imaginary potential obtained from the 
atomic fit, while entering dynamically in the iteration cycles the resulting nuclear density and 
the suppression factor / as defined by Eq.(lO). This dynamical calculation nearly doubled, 
for the light nuclei here considered, the depth of the real part of the phenomenological static 
fT~-nucleus potential of Ref. [3] which is of the 'deep' variety, typically 150-200 MeV deep 
in the static calculation. 

Since there is no preferred way of varying the depth of the real ^"-nucleus potential 
in order to produce different values of binding energies, we used two methods for scanning 
over binding energies. The first one, referred to below as the 'RMF' method, was to scale 
down successively g a x from its initial value g^ K and, once it reached zero, to scale down 
9luK too from its initial value g^ K until the K~ Is state became unbound. As an alternative 

(2) (2) 

method we chose as a starting point values g aK and g^, K obtained from fits to kaonic atom 
data where the RMF-based real potential was joined at large radii by a 'tp' expression, 
using for t the chiral KN amplitudes of Ramos and Oset [15]. This starting potential was 
of the 'shallow' variety [5], about 55 MeV deep in the static calculation. It is gratifying that 
the replacement of the chiral 'tp' expression within the nucleus by the RMF model did not 
change the resulting ^-nuclear potential depth. We then increased the potential depth by 

(2) 

scaling up g C K from its initial value g aK , in order to achieve as deep binding as in the first set 

(2) 

of calculations, while keeping g^K constant at its initial value g uK . This second method will 
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be referred to below as the 'chiral tail' method. Note that in both methods (j = 1,2) the 
starting values {gjjx, 9ax) correspond to potentials that produce in the dynamical calculation 
good fits to the K~ atomic data, although Re V^l have vastly different depths. The good 
fits to the atomic data are inevitably lost once the coupling constants {g a K in the procedure 
outlined above) are allowed to vary in order to scan over a wide range of binding energies 
and the associated widths. 

Figure 1 shows, in its middle and lower parts, calculated widths IV- as function of the 
binding energy Br- for Is states in K l -0 and K 12 G, respectively. Open squares and solid 
circles are for the L and NL versions of the RMF model, respectively, both calculated using 
the 'RMF' method for scanning the binding energies. The crosses are for RMF-L, using the 
'chiral tail' method for scanning the binding energies. It is clearly seen that the widths of 
the K~ nuclear state follow closely the dependence of the suppression factor on the binding 
energy and that, except for binding energies smaller than 50 MeV, the dependence of the 
width on the binding energy follows, for a given nucleus, almost a universal curve. We 
note that the widths calculated in the range B K ~ 100 — 200 MeV assume values 40 ± 5 
MeV, which are considerably larger than what the suppression factor of the upper part 
of the figure would suggest. This is largely related to the dynamical nature of the RMF 
calculation whereby the nuclear density is increased by the polarization effect of the K~ , as 
shown in the next figures. Furthermore, replacing p by p 2 for the density dependence of the 
non-pionic decay modes (2) is estimated to increase the above values of the width by 10 - 
15 MeV. This estimate follows, again, from the increase of nuclear density with B K noticed 
above. Switching on the 7rA decay mode would increase further this estimate by 5-10 MeV 
in the range B K ~ 100 - 180 MeV. We assert that the estimate T K = 40 ± 5 MeV in the 
range Bk ~ 100 — 200 MeV provides a reasonable lower bound on the width expected in any 
realistic calculation. A more detailed systematics is defered to a subsequent regular report. 

Figures 2 and 3 exhibit various nuclear properties for Is states in K l2 G and K 16 0, respec- 
tively. The top and middle parts show the calculated average nuclear density p = \j p 2 dr 
and the nuclear rms radius, respectively, and the lower parts show the Is and lp neutron 
single-particle energies E n . The differences between the linear and non-linear models reflect 
the different nuclear compressibility and the somewhat different nuclear sizes obtained in 
the two models. Again, for K~ binding energy greater than 50 MeV the results are in- 
dependent of the way the binding energy is being scanned. It is interesting to note that 
the increase in the nuclear rms radius of ^- 6 for large values of B K - is the result of the 
reduced binding energy of the lpi/2 state, due to the increased spin-orbit term. Note also 
that as B K - approaches zero we do not recover the values inherent in static calculations for 
the various nuclear entities because the coupling constants g a K and g uK , and the imaginary 
part of the potential, still assume nonzero values. The substantial increase of p, the decrease 
of the nuclear rms radius, and the decrease of the Is neutron single-particle energy, all point 
out to a significant polarization of the nuclear core by the Is K~ . Finally, we note that in 
similar calculations for 40 Ca and 208 Pb the widths of the K~ nuclear state vs. its binding 
energy turned out to be similar in shape to the corresponding results for 12 C and 16 0, but 
the effects on the average density and the rms radius are negligibly small, as expected for 
heavier nuclei. 

Table I shows several examples of the excited-state spectrum of i^~-nuclear states in the 
core of 16 O. The calculated widths of the excited states appear to follow the general trend 
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of widths for Is states as shown in Fig. 1. The polarization of the nucleus as judged by the 
value of its average density p appears to diminish, the higher the K~ state is, in agreement 
with Fig. 3. We note the fairly large spacing between neighboring states: the effective Hcu 
is of order 100 MeV for the first spectrum, decreasing to 80 MeV and to 70 MeV in the 
second and third spectra, respectively, as the Is binding is made lower. Static calculations 
give smaller spacings. Kishimoto et al. [11] have very recently suggested evidence from 
the measured neutron spectrum in the 16 0(K~ ,n) reaction for a peak at Bk ~ 90 MeV 
which they interpreted as the lp K state, probably since a hint for a peak at B K ~ 130 
MeV is also suggested by the same spectrum. The results shown in Table I rule out such 
interpretation which would require Tiuj ~ 40 MeV, considerably less than our dynamical 
calculations produce. If the peak at Bk ~ 90 MeV were a lp state, one should have 
expected the Is state to be deeply bound, at Bk ~ 200 MeV where no signal has been 
observed. On the other hand, if this peak is a Is state, then the relatively shallow lp state 
would be expected too broad to be distinguished in the data. 

IV. CONCLUSIONS 

Binding energies and widths of deeply bound K~ nuclear states in 12 C and 16 O were 
calculated with the aim of establishing values of widths that could be expected for binding 
energies in the range of 100-200 MeV. The method chosen was to couple the K~ dynamically 
to the nucleus within the RMF approach, which is applicable throughout the periodic table, 
but which is disconnected from near-threshold KN phenomenology. Negligible polarization 
effects were found for atomic states, which confirms the optical-potential phenomenology of 
kaonic atoms as a valid starting point for the present study. Substantial polarization of the 
core nucleus was found in these light nuclei for deeply bound K nuclear states. Almost uni- 
versal dependence of antikaon widths on its binding energy was found, for a given nucleus, 
suggesting that the details of how the calculated binding energy is varied over the desired 
range of values is largely immaterial. The widths are mostly determined by the phase-space 
suppression factors on top of the increase provided by the density of the compressed nuclei. 
The present results already provide useful guidance for the interpretation of recent experi- 
mental results [11] by placing a lower limit Tk ~ 35 — 40 MeV on K states in 16 O bound 
in the range B K ~ 100 — 200 MeV. For lighter nuclear targets such as 4 He, where the RMF 
approach becomes unreliable but where nuclear polarization effects are found larger using 
few-body calculational methods [7,8], we anticipate larger widths for K deeply bound states, 
if such states do exist [10]. 

This work was supported in part by the GA AVCR grant IAA1048305 and by the Israel 
Science Foundation grant 131/01. 
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TABLES 



TABLE I. K bound-state spectra in K l -0 calculated for several RMF (NL) Lagrangians 

specified by different coupling-constant ratios a a = g a K/g^K an< ^ a ^ = O^k/o^k- The static 
average density for 16 is p = 0.100 fm~ 3 . 
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FIG. 1. Phase-space suppression factor for the imaginary potential (top), and widths of the Is 
If ^-nuclear state (middle: in 16 0, bottom: in 12 C) as function of the K~ binding energy (see text 
for symbols). 
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FIG. 2. Average nuclear density, nuclear rms radius and neutron single-particle energies for 
j^-C Symbols are as in Fig. 1. 
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